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Abstract 

Scaling relations in four-dimensional simplicial quantum gravity are proposed using 
the concept of the geodesic distance. Based on the analogy of a loop length distribution 
in the two-dimensional case, the scaling relations of the boundary volume distribution 
in four dimensions are discussed in three regions: the strong-coupling phase, the critical 
point and the weak-coupling phase. In each phase a different scaling behavior is found. 
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1 Introduction 



Two-dimensional quantum gravity has generally been regarded as a toy-model theory of four- 
dimensional quantum gravity. According to the Liouville theory, two-dimensional quantum 
gravity can be quantized with a central charge c < 1. The method of dynamical triangula- 
tion |§ has generally been known as a correct discretized model corresponding to the matrix 
model, and has the same critical index as Liouville field theory in the continuum limit. In 
two dimensions many attempts, for example the so-called minbu analysis, j|,|5j the loop length 
distributions, |],[j],[8] the fractal dimensions and an analysis of the complex structures |J reveal 
the many important properties of random surfaces. Especially, we take notice of the concept of 
the loop length distribution. Fortunately, the loop length distribution function has been calcu- 
lated analytically for the case of two-dimensional pure gravity. |7|,§] Our basic strategy in four 
dimensions is that higher-dimensional quantum gravity can be represented as a superposion 
of low-dimensional quantum gravity. [HJ Indeed, a two-dimensional random closed surface can 
be reduced to a summation over the loops along a certain gauge slice regarded as the geodesic 
distance corresponding to time in our case. One of the aims of this article is to investigate 
four-dimensional Euclidean spacetime structures using the geodesic distance (precisely, the ex- 
tended concept of the loop length distribution in two-dimensional simplicial quantum gravity). 
Based on an analogy of the loop length distribution in the two-dimensional case, the scaling 
relations in the four-dimensional case are discussed for three regions i.e. the strong-coupling 
phase, the critical point and the weak-coupling phase. In the three-dimensional case the scaling 
relations of the boundary-area distributions have already been argued using the analogy of the 



loop length distribution, resulting in the establishment of scaling properties. [|Tl|] 

In order to discuss the scaling relations of the boundary volume, we assume that the scaling 

variable (x) has the form V/D a , where V, D and a denote each boundary volume, the geodesic 

distance and the scaling parameter, respectively. 

The rest of the article is organized as follows: section 2 briefly reviews the definition of the 

model. Section 3 considers the new scaling relations of the boundary volume obtained in four 

dimensions, and the following subsections look at these in detail for the three regions. Finally, 

section 4 contains a summary of the results and discussions. 



2 The model 

It is still not known how to give a constructive definition of four-dimensional quantum gravity. 
We naturally must mention the Euclidean Einstein-Hilbert action, as follows: 

Seh = J d 4 x^ (a - ifl) , (1) 

where A is the cosmological constant and G is Newton's constant of gravity. We use the lattice 
action of the four-dimensional model with the topology S A , corresponding to the above action, 
as follows: 

S(k 2 ,k 4 ) = -k 2 N 2 + K A N A 

= -§^+(A , + §«»- 1 (i))^ (2) 

where iVj denotes the total number of z-simplexes and A' = cA; c is the unit volume and 
cos _1 (|) is the angle between two tetrahedra. The coupling k 2 is proportional to the inverse 
bare Newton's constant, and the coupling K4 corresponds to a lattice cosmological constant. 
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Figure 1: Fractal dimensions (df) versus the geodesic distance ploted for the strong phase, 
critical point and weak phase with = 32K. 

For the dynamical triangulation model of four-dimensional quantum gravity we consider a 
partition function of the form 

Z(k 2 ,k,) = J2 e~ s ^\ (3) 

T(S 4 ) 

We sum over all simplicial triangulations, T(S 4 ), on a four-dimensional sphere. Here, we fix 
the topology as S 4 . In practice, we must add a small correction term, |T2[ AS, to the lattice 
action in order to suppress volume fluctuations. The correction term is denoted by 

AS = 5^(N A -Ni tar9et) ) 2 , (4) 
where N^ 9 ^ is the target value of four-simplexes, and 5 = 0.0005 is used. 



3 Numerical simulation and results 

One of the interesting observables in this model is the fractal dimension (df) at large scale. A 
common way to define the fractal dimension is based on studying the behavior of the volume 
within a geodesic distance D, (V(D)). Usually, df is calculated as the following quantity: 

_ d\nV(D) 



which would be a constant if V(D) goes like a power of D. Fig.|I| denotes that the fractal 
dimensions are plotted versus the geodesic distances with N4 = 32K. We have to inquire, 




Figure 2: Number of boundaries versus the geodesic distance for the strong phase, critical point 
and weak phase with double-log scales. 

to some extent, into the internal geometrical view points of the four- dimensional dynamically 
triangulated manifolds. 

We define Nb(D) as the number of boundaries at a geodesic distance D from a reference 
four-simplex in the dynamically triangulated manifolds. Fig.|2| shows the distributions of Nj,(D) 
for the typical three coupling strengths with N4 = 16K. In the strong-coupling limit {n 2 = 0, 
this corresponds to G — > oo), the only boundary that is identified as the mother universe exists 
within almost all of the geodesic distances; also the creation of the branches is highly suppressed, 
which means that the mother universe is a dominant structure. The suppression of branches 
shows one of the characteristic properties of a "crumpled manifold", which is similar to the 
case observed in the two-dimensional manifold. On the other hand, in the weak-coupling phase 
(for example, we chose k 2 = 2.0, (see Fig.|2|)), we see the growth of the branches ("elongated 
manifold") until D ~ 60 for N4 = 16K. From Fig. 1,2 there is no doubt that the dynamically 
triangulated manifold becomes crumple in the strong-coupling phase and a branched polymer 
in the weak-coupling phase. |L3| , [14j 

We now devote a little more space to explaining the loop length distributions in two- 
dimensional dynamically triangulated surfaces. We suppose a disk which is covered within 
a geodesic distance of D from the arbitrary triangle. Because of the branching of the surface, a 
disk is not always simply-connected, and there usually appear many boundaries consisting of S l 
loops in this disk. The loop length distribution function (p(L, D)) is measured by counting the 
number of boundary loops with the loop length (L) within a geodesic distance (D). In ref.7) 
the loop length distribution function (p(L, D)) has been calculated as a function of the scaling 




Figure 3: Loop length distribution function (p(x,D)), consisting of the mother loop and baby 
loops distribution function, (l)mother loop distribution function: yi^e" 1 (thick line) and 
(2)baby loop distribution function: (x~ 5 / 2 + |x _3//2 )e _:r (dotted line). Both are plotted with 
double-log scales. 



variable, x = L/D 2 , in the continuum limit for the case of pure gravity in two dimensions, 

This distribution function consists of two different types of distributions. The first two terms of 
eq.(§p, (x~ 5 ^ 2 + \x~ z l 2 )e~ x , represent the so-called baby loops. The baby loops originate from 
the quantum fluctuation of the surface. 

The last term of eq.(||), (yx 1//2 )e _x , represents the mother loop. The precise definition of 
the mother loop is the boundary of the largest uncovered regions for each geodesic distance. 
One of the notable features of the loop length distribution is that the mother loop and the 
baby loops distribute with the same scaling variable, x = L/D 2 . There are also other things 
to note. The distribution function of the baby loops gives the divergences in the calculations 
of the number of boundaries at distance D, M°(D) = J °° dL p(L,D), and the total length of 
boundaries at distance D, M 1 (D) = J °° dL L p(L,D). These divergences indicate that M°(D) 
and M 1 (D) depend on the lattice cut-off in the model, which offers a key to understanding the 
universality of the distributions. 

For future convenience, we show the mother and baby loop distributions in Fig|| Similar 
distributions also appear in four dimensions! Numerically, these distribution functions were 
investigated in ref.6); excellent agreement between the numerical results and the analytical 
approaches has been obtained. In the sense that two distributions of the sections (or boundaries) 
of the surface at different distances are exactly the same as each other after a proper rescaling of 
the loop length, we can safely state that the dynamically triangulated surfaces become fractal. 
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Figure 4: Boundary volume distribution in the strong-coupling limit with N± = 32K with 
double-log scales. 

That is to say, p satisfies relation (0) under rescaling, D — > D' = VXD and L — > LI = XL, 

p(x,D)=p(L,D) = \- 1 p(L',D'). (7) 

We must now return to the four- dimensional case. The sections (or boundaries) of a four- 
dimensional dynamically triangulated manifold are closed three-manifolds. The boundary vol- 
ume distribution (p(x, D)) is defined based on an analogy of the loop length distribution in two 
dimensions, and assuming that the scaling variable(x) has the form V/D a , where D denotes 
the geodesic distance in the four- dimensional dynamically triangulated manifold. 



3.1 The strong-coupling phase 

Fig|4| shows the boundary volume distributions, p(x) x D 4 ' 5 , with x = V/D A ' 5 as the scaling 
variable in the strong-coupling limit, K2 = 0, while the fractal dimension (d/) reaches 5.5 with 
A^4 = 32K. The boundaries of the mother universe show a good scaling property with x = 
V/D d f~ 1 as the scaling variable, and have a Gaussian distribution. We can recognize that in 
the strong-coupling phase the scaling parameter (a) of the mother universe satisfies the relation 
df — a + 1, and that the manifold resembles a dj-sphere (S f), [plj where df denotes the fractal 
dimension, which increases with the volume [TJj](see Fig.^). 



On the contrary, we have no evidence concerning the scaling properties of the volume of the 
baby boundaries. To be precise, it is impossible to scale the baby volumes in the same manner 
as the mother volume. We find a rapid increase in the distribution of the baby volumes when 
x — > 0. This rapid increase indicates the existence of the "Planck scale" in the model. 
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Figure 5: Fractal dimensions plotted versus the geodesic distance in the strong-coupling limit, 
k 2 = with N 4 = 8K, 16K, 32K. 



3.2 The critical point 

We have long believed that the phase transition in four dimensions is continuous. According to 



recent reports, |15j the phase transition of four-dimensional simplicial quantum gravity seems 
to be first order. Indeed, we have observed the double-peak structure of a histogram of N in 
our measurements at a volume of N± = 32K and k 2 = 1.2581 (see Fig.^), which is a sign of a 
first-order transition, contrary to what is generally believed. We must draw attention to the 
double-peak histogram structure on the critical point. We thus measure the boundary volume 
distribution on both peaks, and obtain the more clear signal of the mother universe on the peak 
which is close to the strong-coupling phase. Fig.[7| shows the boundary volume distribution on 
two peaks. We call the peak close to the strong-coupling phase the "down peak" , and the other 
peak close to the weak-coupling phase the "up peak" . 

The distributions of the boundary volume obtained from the "up peak" are similar to 
the distributions obtained in the weak-coupling phase (see next subsection). That is to say, 
the structure of the mother volume disappears (see Fig|||7|). There is fairly good general 
agreement that the dynamically triangulated manifold behaves like a branched polymer in the 
weak-coupling phase. We may therefore say that the boundary volume distributions of the "up 
peak" are non- universal. 

The data obtained in the "down peak" are shown in Fig.|8|for various geodesic distances(D). 
We observe two different scaling properties: one is for the mother boundary, and the other is 
for the baby boundaries. What we have learned in two-dimensional quantum gravity is that 
the distribution of the mother boundary is universal, and that the baby boundaries are non- 
universal. We should thus notice the mother part of the distribution. Fig.RI shows the volume of 
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Figure 6: Histogram of Vq at the critical point k 2 = 1.2581 with N4 = 32K. 
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Figure 7: Boundary volume distributions of the two peaks with double-log scales, (a) denotes 
the boundary volume distributions of the "down peak" and (b) denotes the boundary volume 
distributions of the "up peak". In the case of (b), the mother volume distribution almost 
disappears. 
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Figure 8: Boundary volume distributions with scaling variable x = at the critical point 
k 2 = 1.2581 with double- log scales. 
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Figure 9: Volume of the mother boundary versus the geodesic distance with double-log scales 
at the critical point k,2 = 1.2581 and 1.2770 with N± = 32K and 64K, respectively. For the 
case of Ni = 32K we estimate the scaling parameter a = 2.4(1) for the range 7 < D < 10. 
For the case of A^ 4 = 6AK we also estimate the scaling parameter, a = 3.0(1), for the range 
7 < D < 12. 



the mother boundary versus the geodesic distance at the critical points K2 = 1.2581 and 1.2770 
with N4 = 32K and 6AK, respectively. In the case of N4 = 32K, we estimate the scaling 
parameter a = 2.4(1)0 for the range 7 < D < 10. In the case of N4 = 6AK we also estimate 
the scaling parameter a = 3.0(1) for the range 7 < D < 12. From Fig.§ we can also roughly 
estimate the scaling variable of the baby boundaries, x = VjD Aa . We cannot measure the 
scaling parameter for the baby boundaries with the same accuracy for the mother boundary. 

In order to discuss the universality of the scaling relations, we assume the distribution 
function in terms of a scaling parameter, x = V/D a , as 

p( x ,D) = ^x a -e^, (8) 

where a\ and 02 are some constants. We obtain the following expectation value of the boundary 
three-dimensional volume appearing at distance D: 

< v®(D) >= - dVVp(x,D), (9) 

IV Jvn 



IV 

where v denotes the UV cut-off of the boundary three-dimensional volume and iV = dV p(x, D) 
is a normalization factor. 

If ai > —1, the right-hand side of eq.(pl) converges with v — > 0, as follows: 



< WHO) >= r '; i + 2 1 ) . g' 

a 2 l (ai + 1) 



(10) 



and gives a finite fractal dimension, df — a + 1. If ax < —1, the right-hand side of eq.@ 
depends on the cut-off t>o, and cannot give a finite value when t> — > 0. In this region, ai < —1, 
we can estimate the dominate part of this integration near to v — > for four regions, as follows: 

' {a 2 log(^)p D a for a x = -1 



Da (ai+2) for _ 2 < ai < -1 

(11) 

v \og(2l)D-« for a 1 = -2 



_ v D- a for ai < -2. 

It is clear that < V^(D) > goes to zero when v — > whenever Oi < — 1. That is to say, the 
distribution with a\ < — 1 is non- universal. 
We also give fractal dimensions by eq.@, 

d f = a + l for - 2 < ai, (12) 

= a + 1 (with logarithmic correction) for a\ = —2, (13) 

and 

d f = 1- a(ai + 1) for a x < -2. (14) 

We can extract the function of p(x, D) from Fig|| and actually find a x — 0.5 for the mother 
volume. As a result, the mother volume distribution is universal. On the other hand, we find 
ax — —2.0 for the baby volumes. It is entirely fair to say that the distribution of the baby 
volumes is non-universal because it depends on the lattice cut-off (vq). A similar non-universal 
distribution of the boundary volume also appears in the weak-coupling phase. 



* A number inside of () is a systematic error from choices of the range of D. The determination of the range 
in which the mother volume expects a power law behavior of D is more or less ambiguous. 




Figure 10: Boundary volume distribution in the weak-coupling phase, k 2 = 2.0 with N 4 = 32K 
with double-log scales. 

3.3 The weak-coupling phase 

In the weak-coupling phase (we chose k 2 = 2.0), as mentioned before, the dynamically triangu- 
lated manifold resembles an elongated branched polymer; thus, we cannot observe the mother 



universe at all. Fig.^Oj shows the boundary volume distributions with x = V/D a.s a scaling 



variable. Moreover, from this figure we can safely state that the p(x) scales as 

p(x) x D oc x- 2 - e~ x , (15) 

and that the distribution is non-universal. 

For all measurements of the boundary volume distributions we prepared about 100 inde- 
pendent configurations on the average and 100 four-simplexes randomly in each configuration 
as the origin of the geodesic distance. 



4 Summary and discussions 

We investigated four-dimensional spacetime structures generated by dynamical triangulation 
using the concept of the geodesic distance. In an analogy of the loop length distribution in 
two-dimensional case, the scaling relations in four dimensions are discussed for the three phases 
i.e. the strong-coupling phase, the critical point and the weak-coupling phase. The boundary 
volume distributions (p(x, D)) give some basic scaling properties on the ensemble of Euclidean 
space-times described by the partition function eq.@. In the case of two dimensions we know 
that the loop length distribution shows excellent agreement with an analytical prediction, |7|,^] 
and that both of the baby loops and the mother loop show scalings with the same parameter 



(x = L/D 2 ). However, since the loop length distribution of the baby boundary depends on the 
lattice cut-off, we can recognize that it is non- universal. 

The question is the boundary volume distributions in the four dimensions. In the strong- 
coupling limit (k,2 = 0) we find that the mother part of the boundary volume distribution 
(p(x, D)) scales with x — V/ D d f~ l as a scaling variable. In this phase there is general agreement 
that the four- dimensional dynamically triangulated manifold seems to be a df-sphere(S df ). 
What is important for the constructive definition of 4D simplicial quantum gravity is the fact 
that the boundary volume distribution of the mother universe in the strong phase gradually 
changes into that at the critical point. The fluctuations of the spacetime growth with K2 — > 
and at the same time the distribution of the baby boundaries comes to scaling (see Figs|| and 

ID- 

At the critical point in the case of four dimensions we have obtained similar boundary 
volume distributions to the two-dimensional case. However, two different scaling relations are 
found: one is for the mother boundary [x = V/D 2A and a\ ~ 0.5) and the other is for the 
baby boundaries (x = V/D 4 ' and a\ ~ —2). For the reasons mentioned above, the boundary 
volume distribution of the mother universe seems to be universal (i.e. there is no dependence 
on the lattice cut-off) and that of the baby universes seems to be non-universal (i.e. there is 
a dependence on the lattice cut-off). We should notice that there exists a considerable size 
dependence of the scaling variable. In the case of A/4 = 32/T and 64AT we obtain x = V/ D 2A ^ 
and V/D 3 '°^\ respectively. As of now, these are all of the data in terms of the scaling variable. 
There is still room for further investigation. 

In the weak-coupling phase we obtained elongated manifolds, in other words, the branched 
polymers. In this phase no mother universe exists and the boundary volume distribution of 
the baby universes shows a scaling relation with scaling parameters a ~ 1.0 and a\ ~ 2.0 (see 
Fig.|T0). However, we recognize that the distribution of the baby boundaries is non- universal. 



We think that the cut-off dependence of the distribution of the baby boundaries in all the 
ranges of K2 is reflected in the existence of the Planck scale in this model. In ref.13), we see that 
the effective curvature increases steeply (diverges) when D —>■ in all phases of k^. 0.8(strong 
phase) < k 2 < 1.5 (weak phase) and the authors suggest the existence of a 'plankian regime'; 
D ~ 7.5 on grounds of the divergence. Our data of the baby boundaries are consistent with their 
view; also, the data of the mother and baby boundaries clear the origin of the scaling relations, 
which has been argued in ref.13), 14). It seems reasonable to suppose that the lattice model 
described by eq.(|3|) corresponds to an effective theory which shows a rapid phase transition 
as well as three different types of scaling relations in each phase. [P~3"|,[r5| In ref.16) we can 
find a discussion of another model which may have a continuous transition. For the present, 
it remains an unsettled question as to how to construct a consistent lattice model including 
"real" four-dimensional gravity in the continuum limit. 

However, our numerical results are expected to be a first step to research the universal 
scaling relations in two-, three- and four- dimensional simplicial quantum gravity. 
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